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We investigate transition form factors of B meson decays into a scalar glueball in the light-cone 
formalism. Compared with form factors of B to ordinary scalar mesons, the B-to-glueball form 
factors have the same power in the expansion of 1/ms. Taking into account the leading twist light- 
cone distribution amplitude, we find that they are numerically smaller than those form factors of 
B to ordinary scalar mesons. Semileptonic B — >■ GW, B — >■ Gl'^l~ and Bs — >■ Gl^l~ decays are 
subsequently investigated. We also analyze the production rates of scalar mesons in semileptonic 
. B decays in the presence of mixing between scalar qq and glueball states. The glueball production 

in Be meson decays is also investigated and the LHCb experiment may discover this channel. The 
sizable branching fraction in Be — ^ {■k^-k~)1~ v, Be — >■ {K^ K~)l~ v or Be — i- {■k^tt' ■k^-k~)1~ v could 
be a clear signal for a scalar glueball state. 
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I. INTRODUCTION 



CN) ■ The existence of the glueball state is permitted by the QCD. Based on the valence approximation, Lattice QCD 
. calculations have predicted that the mass of the lowest-lying scalar glueball is around 1.5-1.8 GeV [l|, 0]. Many 
different candidates have been observed in this mass region, but there is not any solid evidence on the existence of a 
] pure glueball. It is very likely that the glueball mixes with the quark-antiquark state and they together form several 

00 

■ a review on the status. 



physical mesons. On the theoretical side, there are large ambiguities on the mixing mechanism, please see Ref. [3[ for 



Most studies focus on the decay property of the glueball. In fact, the production property of the glueball is an 
alternative way to uncover the mysterious structure of scalar mesons and figure out the gluon component inside. In 
. B meson decays, the color magnetic operator Og^ has a large Wilson coefficient that could produce a number of 
, gluons easily. These gluons in the final state may have the tendency to form a glue state, thus the glueball production 
5^ , in inclusive B decays has attracted some theoretical interests The authors in Ref. 0] have also studied the 

exclusive B — > KG and B — > K*G decays, where G is a scalar glueball. Based on the results in Ref. they have 
neglected gluon recoiled Feynman diagrams and studied the contributions, in which a scalar glueball state is emitted 
and a K or K* meson is recoiled. 

In the present paper, we study the transition form factors of B decays into a scalar glueball and point out another 
interesting mechanism to detect a glueball via the exclusive B decays. These form factors are relevant for productions 
of scalar glueballs in semileptonic B decays. Up to the leading Fock state and leading order in a^, there are three 
different Feynman diagrams shown in Fig.[T] Since the two gluons bounded in a glueball state are already symmetrized 
in the glueball wave function, it is not necessary to consider the crossed diagram. The first diagram is similar with 
the one studied in Ref. 0]. Through the power counting rule, we will show that the first two diagrams have the same 
power with the i3-to-light meson form factors. That will affect numerical results of the transition form factors and 
production rates in semileptonic B decays. 

The ordinary light neutral scalar meson is isospin singlet and/or fiavor SU(3) singlet, while the glueball is flavor 
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FIG. 1: Feynman diagrams of B decays into a scalar glueball G (first row) and an ordinary scalar meson (second row). The 
denotes the possible Lorentz structure arising from the electroweak interactions. 



SU(6) singlet. Therefore it is difficult to distinguish them by the light u, d and s quark coupling. However, the light 
ordinary scalar meson has negligible cc component, while the glueball have the same coupling to cc as that to the uu, 
dd or ss. A clean way to identify a glueball is then through the cc coupling to the glueball. We briefly analyze the 
production in Be meson decays. 

The paper is organized as follows. In section |lTl we give the analysis of the power counting of the _B — > G transition 
form factors. In the perturbative QCD (PQCD) approach 0,0], the form factors arc calculated in the section Hill 
Numerical results and detailed discussions will be presented in the subsequent section. We also briefly analyze the 
discussion in Be meson decays. We conclude in the last section. Functions used in the PQCD approach are collected 
in the appendix. 



II. POWER COUNTING OF TRANSITION FORM FACTORS 

Transition form factors of B decays into scalar mesons are defined by 
{SiPs)\ql^.l5b\BiPB)) = - 
{S{Ps)\qa^,b\B{PB)) = ie^,,o.0Psq^ 



{Pb + M - 2 
2FT{q') 



F,iq')+ "'^ ,' q,Foiq' 

q 



rriB + ms 



{S{Ps)\qa^,j5b\B{PB)) = ^^^[{PB + Ps)^.qu-q^.{PB+PsU, (1) 

niB + ms 

where 5* denotes an ordinary scalar meson or a scalar glueball. q = Pb — Ps and niBims) is the mass of the B (scalar) 
meson. 

In the following, we will work in the rest frame of the B meson and use light-cone coordinates. In the heavy quark 
limit, the mass difference between b quark (rrib) and B meson (m^) is negligible. The mass of a light glueball is small 
compared with the b quark mass, thus in the transition amplitudes we keep them up to the leading order. Since the 
scalar glueball in the final state moves very fast in the large-recoil region, we choose its momentum mainly on the 
plus direction in the light-cone coordinates. The momentum of B meson and the light scalar meson can be denoted 
as 

Pb-^(1,1,0^), P2 = ^(P,0,0^). (2) 
The momentum transfer ia q = Pb — -P2, and there exists the relation p sa 1 — q^ /m^. 
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In the expansion on as , the lowest order Feynman diagrams for form factors of B decays into a scalar glueball are 
depicted in Fig. [T] Up to the leading Fock state, two gluons should be generated to form a glueball. In exclusive B 
decays, these two gluons can be emitted from either the heavy b quark or the light quark. The (8) in Fig. [T] denotes 
the weak vertices. Momentum fractions of the light antiquark in B meson and the upper (lower) gluon depicted in 
Fig [T] arc denoted as xi and X2 {x2 = 1 — X2)^ respectively. 

In the first diagram of Fig. [1] each of the two internal quark propagators is the sum of a collinear momentum 
(gluon) and a soft momentum (light quark). The virtualities of them are of order AgcD^h, where Aqcd is the 
hadronic scale. In the second and the third diagram, one or two light propagators become heavy b quark propagators, 
whose virtualities become instead of AgcD'Tib. Superficially, one may conclude that the power counting for the 
three diagrams obey the relation: F{a) > F{b) > F{c). But in fact, this relation is not exactly correct. As we 
will show in the following, the leading twist light-cone distribution amplitude is constructed in the case that the two 
gluons are transversely polarized. So the structures of the vertices attaching to these two gluons in the first diagram 
have the form: 7^. Apparently, the numerator of the propagators commutes with the transverse Dirac matrix. So, 
the amplitude is proportional to (a;2-^2 — xiIfB){If2 — xilfs) ~ —xilfslfi — xxX2lf2lfB- We can see that neglecting 
the glueball's mass square, there is a small momentum fraction xi in the numerator, which will cancel one of the 
momentum fraction of the denominator. The effective power for one light propagator in the first diagram becomes 
771^, which is the same as a heavy propagator in the second diagram. It implies: F{a) ~ F{b) > F{c). Adopting 
the power counting rule for the B meson wave function and the scaling behaviors for the distribution amplitudes of 
collinear meson, given in Ref. [lo| . we directly obtain 

~ asi^^^^J^) (^^y , (3) 

where the form factor is dominated by the first two diagrams. The B ^ G form factors have the same scaling rule 
with the i?-to-light transition form factor. And we can expect that the gluonic S — > 77 form factors also obey this rule. 
While the light-cone distribution amplitude of a gluonic pseudo-scalar meson is normalized to zero, only the higher 
Gegenbauer moments contribute. The first effective Gegenbauer moment of rj and 77' is very small, so the gluonic 
contribution to B ^ rj form factors is found to be numerically small 0] • 

Since these diagrams are free of the endpoint singularity, both collinear factorization and kx factorization are 
applicable. In the following, we will use the perturbative QCD approach to investigate the transition form factors, 
where the inclusion of the transverse momentum only improves the theoretical precision. 



III. TRANSITION FORM FACTORS IN THE PERTURBATIVE QCD APPROACH 



B meson is a heavy-light system, whose light-cone matrix element can be decomposed as 
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(27r)' 



-e^''^-^OMQ)qa{z)\B{PB)) = 



{pB + ?77-b)75 



(f)B{ki 



V2 



(pB{ki 



(4) 



pa 



where n = (1,0, Ot) and v = (0,1,07-) are light-like unit vectors. There are two Lorentz structures in B meson 
light-cone distribution amplitudes, and they obey the normalization conditions: 



(j)B{ki 



B 



0B(fcl) =0, 



(5) 



with /b as the decay constant of B meson. 

In principle, both the (pBiki) and (pB^ki) contribute in B meson decays. For i?-to-glueball transition, the (j}B{ki) 
can be neglected in the second diagram. The amplitude for the form factors is obtained by evaluating the trace of 
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the wave functions and the scattering kernels. As we will show in the following, the gluon in a glueball is transversely 
polarized at the leading power. So the Lorentz structure 'tfi for the wave function (t>B{ki) commutes with the gluon 
vertex and encounters the propagator [xnlf2 — x\^b)- The momentum of the glueball is parallel with this direction n, 
thus the nonzero term contains a factor x\ and is naturally power suppressed. For the other diagrams, we do not find 
any analytical reason to neglect the distribution amplitude 0b (/ci). The computation depends on the shape of this 
distribution amplitude but it is not well constrained in the PQCD approach at present. In fact, including the ^^(fci) 
term will not improve the quality of the calculation significantly, but introduce one more free parameter. Nevertheless 
in B decays into ordinary qq mesons the contribution of 4>B{ki) is usually neglected, because its contribution is 
numerically small in the PQCD approach 0, [3- We will keep the term with (/)s(fci) in equation ^ in this work 
and leave the term with 4>B{ki) in a future study. In the momentum space the light cone matrix of B meson can be 
expressed as: 



V6 



[fs + m_B)750B(fci). 



(6) 



Usually the hard part is independent of fc+ or/and k , so we integrate one of them out from (j)B{k^ , k , kj^). 
For the light-cone wave functions of the B meson, we use 



^B,apix,b) 



/2Nr 



[PB15 + "ls75]„/3 'f>B{x, b), 



(7) 



where x is the momentum fraction of the light quark in B meson and 6 = |b| , which is the conjugate space coordinate 
of kx. In this paper, we use the following model for b): 



4>b{x, b) = Nbx^{1 - x)^exp 



2^1 



\{^bbf 



(8) 



The normalization factor Nb is determined by normalization condition. 
Decay constant of a scalar glueball state is defined as: 



(9) 



where F^i, is the gluon field strength tensor, fg is determined as = (100 ^ 130) McV [15j. When the two gluons 
are separated in coordinate space, the nonperturbative matrix elements are parameterized in terms of the light cone 
distribution amplitudes (LCDAs). According to the conformal symmetry, the fundamental field with definite twist 
is the component of the gluon field strength tensor 16[. The component F^^ is twist-1, F^^ and F^ are twist-2 
while F~-^ is twist-3. In the following, we only consider the leading twist light-cone distribution amplitudes of the 
glueball state: 



{G{P)\F''^+^'{z-)F^'+•'{Q)\Q) 



dxe 



- "'^^'^(^^K>G(X)] 



(10) 



where n, v are the transverse indices. The coordinate z has been chosen on the light-cone z^ = 0. In Eq. (fTO|) we 
have used the light-cone gauge so that the gauge links between the field strength tensors vanish. The distribution 
amplitude (j){x) can be expanded in terms of the Gegenbauer polynomials: 



0G(a;) = 30x'^{l~xf 



l + ^a„C^/2(2x~l) 



In the ~ gauge, the light cone distribution amplitude is reformulated as: 
{GiP)\A'^^{z-)A''^mO) = -^'dxe"^-^+-^^!^ 



(kcix) 



x{l ~ x) 



(11) 



(12) 
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where only the transverse gluon contributes at the leading twist. 

Before presenting the formulas for the transition form factors, we will briefly comment on the transverse distribution 
of the wave functions. The basic idea of the PQCD approach is that it takes into account the transverse momentum 
of the valence quarks in hadrons which kills the endpoint singularity. The form factor is expressed as a convolution of 
the wave functions and a hard kernel. Rcsummation of the double logarithm due to higher order corrections results 
in the Sudakov factor. Strictly speaking the transverse distribution in all wave functions should be taken into account 
in the PQCD approach. One of the most acceptable candidates would be the exponential wave function or improved 
exponential form like the Gaussian one for the B meson given in Eq. ([8]) . The common feature of these distributions 
is that the contribution from the large b region is exponentially suppressed. Meanwhile the Sudakov factor can also 
suppress the contribution from the large b region as shown in Fig. 2 in the first paper of Ref. Q. This suppression 
effect also depends on the longitudinal momentum. Since the momentum of the quark and antiquark in a light meson 
is large, contributions from the large b region are strongly suppressed by the Sudakov factor. As a consequence the 
role of the transverse distribution in the wave function for the final mesons has already been fulfilled and it can be 
neglected in the PQCD approach. On the contrary, the Sudakov suppression for the B meson is not so manifest that 
the transverse distribution is required. The commonly-used form in the PQCD approach is the one given in Eq. (|8]) 
which will be adopted in this work. 

The form factors in the PQCD approach are obtained by evaluating the three diagrams in Fig. [TJ They can be 
expressed as the convolution over the longitudinal momentum fraction and the transverse space coordinates: 

Fi{q^) = 'i ^'lf '^^^ 'KCFm\ ( dxidx2 [ bidbicj)Bixi,bi) '^'^^^'^^ 

N^~l Jo Jo X2{l-X2) 

] 

62^62 {a^l [p + (1 - P)X2] Ea{ta)ha + pX2Ec{tc)hc} + pX2{l - X2)Eb{h)hb > , (13) 

4—5 -TTCprng / dxidx2 / bidbi(pB[xi,bi) — - 

N^-l Jo Jo X2{l-X2) 

X |y b2db2 {Xip [2- p+{p- 1)X2] Ea{ta)ha + p{2 - p)x2Ec{tc)hc\ 

+p'x2{l-X2)Eb{U)hb] , (14) 

Friq ) = 4 — ttCfTOb / dxidx2 / 6id&i(^s(xi, &i) — 

- 1 Jo Jo 3:^2(1 -2:2) 

b2db2 {xiX2Ea{ta)ha + a;i(l - X2)Ec{tc)hc} - XiEb{tb)ht > , (15) 







where Cp = 4/3 and — 3 are color factors. Ei,hi are hard functions determined through the propagators in 
the three diagrams, which are collected in the appendix. Through fitting the results in the hard-scattering region 
< < lOGeV^, we extrapolate them with the dipole model parametrization 

^^^'^'^ " l-a(gVm|f+6(qVm|)2' ^^^^ 

where i = 1, 0, T and a, b are parameters to be determined in the fitting procedure. 

In the numerical calculation, we adopt ojb ~ (0.40±0.05)GeV and fs = (0.19± 0.02)GeV for B mesons Q. For the 



Bg meson, we use ws, = (0.50± 0.05)GeV and fs^ ~ (0.23± 0.02)GeV [17[. The decay constant of the scalar glueball 



is used as = 0.13GeV [15|, but so far there is no theoretical study on the LCDAs of the glueball. The simplest 
choice is to use the asymptotic form. To roughly estimate the uncertainty from the higher Gegcnbauer moment, we 
will also try the Gegenbauer moment 02 ~ 0.2. The results in table |I] show that the form factors arc sensitive to the 
Gegenbauer moment 02, which will provide relatively important uncertainties to our results. The mass of the glueball 
is taken as mc = 1.5 GeV. In the PQCD approach, the i?-to-glucball transition form factors arc not very sensitive to 
this mass. 
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TABLE I: Distinct contributions to B — t- G form factors at = 0: tlie index (a,b,c) denotes the contribution from the diagram 
(a,b,c) in Fig. [1] respectively. Except the results in the last column where 02 = is used, all other results are obtained with 
a2 = 0.2. 





a 


b 


c 


total 


asymptotic 


B Fi^Fo 


0.042 


0.012 


0.001 


0.055 


0.043 


Ft 


0.035 


-0.011 


0.001 


0.025 


0.017 


Bs^G Fi = Fo 


0.038 


0.011 


0.002 


0.051 


0.039 


Ft 


0.032 


-0.012 


0.001 


0.020 


0.014 



TABLE 11: B ^ G and Bs ^ G form factor with the dipole parametrization. Results in the "asymptotic" ("total") rows 
correspond to the Gegenbauer moment a2 — (02 — 0.2). 





Jb(0) = Fi(0) Ft(0) a(Fi) 6(Fi) a(Fo) b(Fo) a{FT) b^Fr) 


B ^ G asymptotic 
total 


n ri/l 0+0. 008+0. 004 00,7+0.004+0.001 , 99+0.00 n 1 c+O.OO n oc+0.01 n f\a+0 02 , cq+O.OI n cn+O.Ol 
'-'•"^•5-0.007-0.004 ^•"^'-0.004-0.001 ^'^^-0.02 "-^^-0.02 "•""-0.02 "•'-'"-0.03 ^^"^-0.03 "•"^-0.06 

n nt=;t^+o-oii+o.oo5 o O9r+o.oo6+o.oo2 , 9^+0.01 017+0.01 n 01+0.00 A nq+o-oi 1 R9+0.01 0^9+0.02 

"•"'-''-'-0.009-0.004 "•"^'-'-0.005-0.002 ^^^^-0.02 "•^'-0.02 "•"^-0.02 "•"'-'-0.04 ^^"^-0.04 "•"^-0.08 


Bs ^ G asymptotic 
total 


n noQ+0. 006+0. 004 o oi .+0.003+0.001 i 90+O.OI 1 q+0.02 o/|+0.02 o^+0.03 1 a/i+0.02 p. 
0.039_o.oo5-0.003 O.Oi4_o oo3_o i.2U_o oi U.i3_o oi U.84_o oi O.U4_o oo i.D4_o o2 O.D3_o (,2 

n nt^l +0 008+0 005 n 090+0 005+0.002 1 99+0.01 oir+0.02 70+0. 02 n ni +0.03 1 r7+0.02 n r/i+O.QS 
U.Udi_o. 007-0. 004 '-'•"^"-0.004-0.002 ^^^^^-0.01 "•^^'^-0.01 "•'^-0.01 "•'-'^^-0.01 ^^•'^'-0.02 '-'•'^^-0.01 



The distinct contributions from the three diagrams are collected in table ID Among them, the largest contribution 
is from the first diagram, while the third one is smallest. This is in agreement with the power counting rule given 
in Sec. |lTl Although the third diagram can be neglected compared with other potential contributions like the higher 
order corrections and higher power corrections, this diagram is taken into account since it has the same order in as 
with the other diagrams. Compared with the gluonic contributions to the S — > ?/ form factors in Ref. our results 
are much larger. For the pscudoscalar meson 77 and ?/, the gluonic LCDAs are normalized to 0, thus only higher 
Gegenbauer moments contribute. The Gegenbauer moment used for the gluonic content of 77 and 77' is very small 
i?2 = ^ ~ -gg-- Moreover, we have taken into account the second diagram, which also makes the form factors larger. 

In the following, we will comment on the magnitude of i3-to-glueball form factors and compare with other transition 
form factors. In the PQCD approach, one ingredient for the form factor is the matrix element of the giuon fields and 
the quark fields. These matrix elements differ with the LCDAs by some constants which can be explicitly given in 
the definition of the LCDAs. For the glueball state the matrix element associated with the asymptotic twist-2 LCDA 
defined in Eq. ((T2|) can be rewritten as 

1 (j)G{x) _iQx{l~x) 
2{N^ - 1) x{l - x) ~ 2(iV2 - 1) ' 

where we have substituted the asymptotic form for the LCDA (j)Q{x). Similarly, the matrix element defined by the 

TABLE III: B — >■ /o(1370) and Bs — >■ /o(1500) form factors in two scenarios (denoted by SI and S2) with the dipole parametriza- 
tion, where /o(1370) is assigned as fin and /o(1500) is identified as ss. 





Fo(0) Ft(0) a(Fi) b{Fi) a{Fo) b{Fo) a{FT) b^Fr) 


B ^ /o SI 
S2 


n QA+O OS _n QO+O-IO 1 co+0.09 o ro+0.14 o ^0+0.07 _o 91+O.I5 1 cn+O.Oe rQ+0.08 
— U.,3U„o 09 '-'•'5^-0.11 ^^•"'J-0.05 "•"'J-0.08 "•"-'-0.02 "•^^-O.OS ^^•"'-'-0.04 "•"'-'-0.05 

0.63l«:f, 0.76l°:f, imtrol 0.53^-^ 0.70l°:?? -O.Ut'^f, 1^63««^ 0.57t°«^ 


Bs fo SI 
S2 


-0.26lH^8 -0^34t«:i° 0.72trot -0-20tVi'o 1^6llS:J^3 OAStUl 1^60t^:^! 0.48««^ 

f.o+0.20 89+0-30 n fi^+0 04 o 99+O.O7 17C+O.13 71 +0.20 1 71 +0.04 ^(.+0.06 
'-'•'-'"-0.12 U.OZ_o 15 U.DO_o.io — U.^^_o.02 ^•'D-0.08 '^•'J^-O.OS ^-'^-0.07 "•OD-0.10 
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quark fields for a pseudoscalar meosn is given by 

4>p{x) 6x{l — x) 
~4JV7 " 4Nr_ ' 



where the asymptotic form for the twist-2 LCDA 4>p{x) is used: (j)p{x) = 6x{l — x) [16j. The prefactor for the glueball 
is almost 4 times larger than that for a pseudoscalar meson. Despite it, the i?-to-glueball form factors are smaller 
than the B — >■ tt form factors. Compared with our previous studies [3 [3 on the transition form factors of B meson 
decays into ordinary scalar mesons in Tab. lIIH the B ^ G form factors are smaller by a factor of 6— 10. One reason is 
that the decay constant for glueball is smaller than the scalar decay constant for the quark content. For example, the 



scalar decay constant for ao(1450) is —0.28 ±0.03 or 0.46 ±0.05 [20|. There exist some other reasons. In the transition 
form factors for the quark content, the small momentum fraction xi in the numerator of the first diagram in Fig. [1] 
is replaced by the factor accompanied with the twist-3 LCDA: 2rs = Although all of these two terms have the 

same power, the transition for the quark content is much larger than that for the gluon component since 2rs 3> (xi) 
((xi) is the typical value of xi). For the second diagram, the longitudinal momentum of the quark propagators are 
given as: X2X2X1, while the corresponding one that appears in the quark transition form factor is x\xi. In the latter 
case, the region with small values of X2 will give relatively large contributions although the endpoint singularities are 
removed by the transverse momentum and the threshold resummation. In the former transition, the form factor is 
directly reduced by the propagators and the contributions from the small momentum fraction region are small. The 
third diagram is power suppressed as we have discussed. Thus the total form factors are smaller than those for quark 
contents. 

Results for the B ^ G form factors, together with the (/^-dependent parameters in Eq. (fT6| . are collected in Tab.llll 
The uncertainties are from the input parameters: (i) the B meson decay constant and the shape parameter (wb,^^^) 
in the wave functions; (ii) the factorization scale from 0.75i to 1.25t (not changing the transverse part l/bi): 

max [{^.Tb^/pximB, l/^^i, I/&2] <ta< max [1.25^pxims, l/^i, I/&2] 7 
max [0.75y^px2"ib, y/ px2XimB, I/&1] < tb < max [l.25y^px2mB, y/ px2XimB, 7 
max [0.75-y/pTOB, I/61, 1/62] <tc< max [1.2b^mB, l/f*!, 1/^2] 7 

and the hadronic scale Aqcd = (0.25 ± 0.05) GeV. The g^-dependence of the B — > G form factors, together with the 
form factors for ordinary scalar mesons, are given in Fig. [51 The dashed (black) and solid (red) lines denote the form 
factors for ordinary mesons in scenario I and scenario II, while the dash-dotted (blue) lines denote the _B — > G form 
factors. 

Although form factors of B decays into scalar glueballs are smaller than those for the quark content, one can not 
neglect the gluon content in the case of mixing for scalar mesons. There is a nontrivial factor l/\/2 in form factors for 
the quark content in B meson decays. This factor will make the glue component in a scalar meson more important. 
Compared with the ordinary B ^ S form factors, the PQCD calculation of the B ^ G form factors is expected to be 
more reliable. In the ordinary B ^ S form factors, the perturbative hard-scattering diagrams contain the endpoint 
singularity in the collinear factorization. Although the inclusion of the transverse momentum can resolve this problem, 
the results would still be sensitive to the treatment of the endpoint region. In the PQCD approach resummation of 
the double logarithms results in the Sudakov factor which will suppress contributions from the nonperturbative region 
(large b region). This will improve the convergence of the perturbation theory. These Sudakov effects in — > tt 



transition have also been investigated in Ref. [2l| with a different conclusion (see Ref. [221 for the response). On the 
contrary, the situation is different in the i?-to-glueball transition which is free from the endpoint singularity. The 
results are more stable than the B ^ S form factors. 
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FIG. 2: The (/^-dependence of B S form factors and \Fo{q^)\, where the dashed (black) and sohd (red)hnes denote 

the form factors for ordinary mesons in scenario I and scenario II, while the dash-dotted (blue) lines denote the B ^ G form 
factors. 



IV. PHENOMENOLOGICAL APPLICATIONS 



A. Semileptonic B — s- Glf and B — s- Gl'^l decays 



„ 10 

Hcff = ^^VtbV;^ (18) 



Integrating; out the highly off-shell intermediate degrees of freedom, the effective Hamiltonian for b — > ulh'i transition 
is given by [23| 

n,s{b ^ uWi) = ^Kf,S7M(l - 75)W7''(1 - IrM- (17) 
while the effective Hamiltonian responsible for the b Dl^l^ [D — d, s) transition is 

Gp 

where Vut, Vtb and Vto {D = d, s) are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. Ci{fi) are the Wilson 
coefficients and the local operators Oi{^) are given by 

Oi = {DaCa)v-Aicf}bp)v-A, O2 = (riQC/3)y _ A (c/3&a )y- A , 

O3 = {Daba)v-AYiqi3qi3)v-A, Oj = {D abp)v -A'^Xq^qc^V -A, 
9 <? 

O5 = {Daba)v-A'Y{qi3qi3)v+A, Op = (Z)a&ff )v-A y^(gffga)v+A, 

f?5.-(l!,,)5^,. + ^- 
^(;V)(^7'^(l-75)&), 0,, = ^ 



O7 = ^^Da^^'^il + ^,)bF,, + -ffDa^'il - l,)bF,,, 

O9 = ^MiDri^ - 75)&), Oio = ^(;7;.750(^^7^(1 - 75)6), (19) 



where (gig2)y-A(93'74)y-A = (9i7''(l - 75)92)(g37A<(l ^ 7)94), and (gig2)y-A(g'3'74)y+A = (q'i7''(1 - 75)92)(g37M(l + 
7)94). a,/3 are the color indices for the quark field, aom = 1/137 is the fine structure constant. The amplitude for 
b — > sl^l^ transition can be decomposed as 

A{b ^ si+n = ^^vttv*l^cf{q')[s^,{i - i.miri] + cMsi,{i - iMn^i] 

-2mfaCf J(l + 75)6] [^7'^^] - 2msCf[sia^,^{l - 75)6] [1^1]^, (20) 
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TABLE IV: The values of Wilson coefficients Ci(jnb) in the leading logarithmic approximation, with mw = 80.4GeV, ^ 
mi,,poie 0- 





Ca 


Cs 


Ca 


Cs 


Cs 




Cio 


1T07 


-0.248 


-0.011 


-0.026 


-0.007 


-0.031 


-0.313 4.344 


-4.669 



with mb as the b quark mass in the MS scheme and we will use rrib = 4.8 GeV. Cf^ = C7 — C5/3 — Ce and Cg^ 
contains both the long-distance and short-distance contributions, which is given by 



(21) 



with s = /m\. Import represents the perturbative contributions, and Ild is the long-distance part. Since the long- 
distance contribution can be separated on the experimental side, we will neglect it. The short-distance corrections 
Fpcrt is given by 
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Fpert(s) - Mm„s)Co--/i(l,J)(4C3 + 4C4 + 3C5-f Ce) 

-i;i(0, 5)(C3 + 3C4) + ^(3C3 + C4 + 3C5 + Ce), 



(22) 



25|. The relevant Wilson 



with Co = Ci + 3C2 + 3C3 + C4 + 3C5 + Cq and = rric/mb with = 1-27 GeV 
coefficients, listed in table ITVl are given up to the leading logarithmic accuracy [2^. 

Lepton decay amplitudes can be directly calculated using the perturbation theory. The unknown amplitude is the 
matrix elements of quark operators which have been parameterized as the form factors. The partial decay width is 
given by 



dT{B Glv) 



dq^ 



V\Gl\Vub? fq^ 
3847r3m|(72 



,2\ 2 



{mt + 2q')XFt{q') + 3mf (m^ - misYF^{q') 



(23) 



where A = {m?g 



™g) ^ ^^m^q and mi is the lepton mass. Integrating out the momentum transfer square g , 



we obtain branching ratios of _B — > Glv (l = e, fi) and B — > GtVt 

B{B^Glv) = (0.18 
B{B Grvr) = (0.08 



, -1-0.05-1-0. 02-|-0.01i\ „ in- 
-0.03-0. 02-0. 01J-' ^ 



-1-0. 08-1-0. 03-1-0. 02 rr, -i -i H 
-0.06-0.03-0.02 [^^-^ ^- 

-1-0. 03-1-0.01-1-0. 01 [n nc;+0-02-l-0.01-l-0.00i>, ^ 2^0-5 



O.Oi-l-O.Oi rr> r.r-|-0.0i!-|-0.0i-|-0.0 

0.03-0.01-0.01 ['-'■"'3-0. 02-0. 01-0. OOJ 



(24) 
(25) 



where results in the square brackets and out of the square brackets are evaluated using the asymptotic form and the 
form with 02 ~ 0.2, respectively. In the evaluation of semileptonic B decays we have adopted the dipole form for the 
form factors given in Eq. (|16p . where the parameters are obtained through fitting the large recoil region. Uncertainties 
are from three different kinds of inputs: (i) the B meson decay constant and the shape parameter (ijJb^<^b,)-^ (ii) the 
hadronic scale Aqcd, the factorization scale i, (iii) the CKM matrix element \Vub\ 
with the recently measured result of _B — > i]li> decay 27 1 



(3.5ltJ^:ig)xl0-3 [261. Compared 



BiB^^rjl-i?) = (3.1 ±0.6 ±0.8) X 10" 



(26) 



we can see that B — > GW{1 ~ e, /i) have smaller branching ratios by one order. The main decay channel of a scalar 
glueball could be tttt or KK while the rj meson is reconstructed by three pious or two photons. Thus a scalar glueball 
is easier to detect than the isoscalar meson 77, and the B — > GW decays may be observed in the future. 
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The partial decay width of semileptonic Bd.s — ^ Gl'^l decays is given as 



dTjBd,, ^ Gl+l- 



X IVtbVtD] 



1024m|j7r 



iriB + mp 



2 2 y (ms + mp) 



4Xmf 



(ms + Top) 



+ ^|Cio(TO|-77i|,)Fo(g2)| 



(27) 



where D = d,s for i?'' Gl^l and — > Gl^l decays, to; denotes the lepton's mass and Gtl 
^^^^^Gj^ ^ 0. Branching fractions of _B ^ Gl^l^ and Bg Gl^l^ decays are predicted as 



B{B 

B{B- 

B{Bs 
B{Bs- 



^ Gi+r] 

Gt+t-] 

^ Gi+r] 

Gt+t-] 



— fn 07+0.12+0.05 [n 1 0+0. 07+0. 03iN „ -in-o 

— ^U.Z(_Q Qg_Q Q4[U.iO_Q Q5_Q Q3JJ X iU , 

= (o.5l™:i[o.3tEJ:liri])xio-", 
= (o.6t™;l[o.n.ii"o:l])xio-«, 

+0.07+0.04 jQ ^2 



(0-1^-0.05-0.03 



+0.04+0.02] 
O.O3-O.O2J 



X 10" 



(28) 
(29) 
(30) 
(31) 



Uncertainties are from two different kinds of inputs in the PQCD approach: (i) the B meson decay constant and 
the shape parameter (wsjOJbJ, (ii) the hadronic scale Aqcd, the factorization scale t. The CKM matrix element 



Vtd = {8.59toil) X 10-3, Vts = -iA0A2tolj) x lO'^ and Va = 0.999146t[!:[!!][J!]^^ [261 wih not give large uncertainties. 
The B — > Gl~^l~ has tiny branching fractions, which can not be observed on the present B factories, especially the 
B Gt^t^ with a very small phase space. But the Bg Gl^l^ decay channel may be observed on the future 
experiment, since it has a sizable branching fraction. 

The physical scalar meson is probably a mixture of glueball and ordinary states, so the predicted branching ratios 
in eq. p4|25|28ll3T|) will be smaller for physical scalar meson, according to the gluon component of each scalar meson. 
The experimental accessability will still be large, at least for the scalar mesons with large glue components. 



B. Mixing between scalar mesons 

There is not any solid experimental evidence for a pure glueball state up to now. Lattice QCD predicted the mass of 
a scalar glueball ground state around 1.5-1.8 GeV . At present most of the lattice studies on hadronic correlators 
are in the quenched approximation by neglecting the fluctuations of the quarks. Due to the lack of dynamical quarks, 
the binding of quark- ant iquark systems must be attributed to the nonperturbative properties of gluons, the unique 
dynamical degree of freedom in the Lattice QCD. Secondly the simulations are based on the computation of the 
matrix element of the gluon operators, while the glueball states G are obtained by smeared gluonic operators since 
there is no physical glueball state. Nevertheless, despite these ambiguities the Lattice QCD simulations give us a hint 
that one scalar glueball state might exist around this mass region. It is very likely that the glueball state mixes with 
the ordinary quark-antiquark state to form several physical mesons. 

In the literature, three scalar mesons /o(1370), /o(1500) and /o(1710) are usually considered as the potential 
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candidates. The mixing matrix can be generally set as 

/o(1710)\ / ai 02 as 

/o(1500) = 61 ^2 ^3 1 I I . (32) 
/o(1370)/ \ci C2 C3 

The unitary condition implies that the mixing matrix has only three independent real numbers. A non-zero ai, 61 or 
ci would be a clear evidence for the existence of a glueball. The semileptonic B JqIv decays receive contributions 
from the fin component and glue component but without ss component (at least negligible), while the semileptonic 
Bs — > fol'^l^ channel only receive contributions from the ss and glue component but without nn component. Thus 
the experimental measurements for the two channels can give effective constraints to the mixing parameters. For the 
three kinds of flavor singlet scalar mesons, we have altogether 6 experiments, but only three real parameters in eq. (|32p 
to be fixed. Since the branching fraction of Bs fol^l^ is expected to have the order of 10^® or even smaller, one 
needs to accumulate a large number of B decay events, which could be achieved on the future experiments such as 
the Super B factory. 

Taking into account the current available experimental data, the mixing mechanism is still not fixed, for example 
there are two quite different mixing matrix determined by two group of people. Because the decay width of /o(1500) is 
not compatible with the ordinary qq state, Amsler and Close claimed that /o(1500) is primarily a scalar glueball j28| . 



In the subsequent studies 29| , they extracted the mixing matrix through fitting the data of two-body decays of scalar 
mesons 

/ /o(1710) \ / 0.36 0.93 0.09 \ f G \ 

/o(1500) = -0.84 0.35 -0.41 ss . (33) 
V .fo(1370) / V 0.40 -0.07 -0.91 J \nn J 

Based on the SU(3) assumption for scalar mesons and the quenched Lattice QCD results Q, Cheng et al. 30 1 reanalyze 
all existing experimental data and fit the mixing coefficient as 

/ /o(1710) \ / 0.93 0.17 0.32 \ ( G \ 

/o(1500) = -0.03 0.84 -0.54 ss • (34) 
V .fo(1370) / V -0.36 0.52 0.78 ) \ nn ) 

Here, the /o(1710) tends to be a primary glueball. This is very different from the first matrix of mixing coefficients 
in p3p . The scalar meson production rates in B meson decays can be used to distinguish these assignments, starting 
with the B ^ S form factors collected in Tab. HI] and Tab. IIIII 

• In scenario I, if we use the mixing coefficients in Eq. (|33p . the production rates of /o(1710) and /o(1500) in 
B decays are much smaller than that of /o(1370) but they have large and comparable production rates in Bg 
decays; if we use the mixing coefficients in Eq. /o(1710) has small production rates in both B and Bg 
decays but the other two mesons have large and comparable production rates in B and Bs decays. 

• In scenario II, if we use the mixing coefficients in Eq. ([33|) . production rates of /o(1370) and /o(1500) in B 
decays are large and comparable, and /o(1710) is copious in Bs decays; if we use the mixing coefficients in 
Eq. (p4l) . three mesons have similar production rates in B decays, and /o(1500) is more copious than /o(1370) 
and /o(1710) in Bs decays. 

Based on our predictions on form factors in this work and in our previous studies [isi [l^ , these differences in B and 
Bs decays are helpful to distinguish the two mixing mechanisms. Once the branching fraction of these decays were 
measured, our calculation can be used to constrain the mixing angle. One more practical method in future would 
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be to use the x^-fit method and take the theoretical and experimental uncertainties into account. This requires the 
future experimental studies. 

The main decay channel of a scalar meson would be G ^ PP, and thus the scalar meson could be reconstructed 
from TT+TT^ or K^K^ . The scalar meson may also decay into a pair of vectors and thus it could be reconstructed as 4 
pious. For the three candidates of glueball states, the reconstruction method is different. If /o(1370) is the dominant 
glueball state, the best candidate is from 4-pion state, since its main decay mode is /o(1370) — > pp [2^. For /o(1500), 
the dominant mode is /o(1500) 7r7r(34.9 ± 2.3%) and /o(1500) 47r(49.5 ± 3.3%). As for /o(1710), the dominant 
decay mode is /o(1710) ^ KK. 



C. Uncertainty analysis 

Besides the uncertainties that we have already given in the above, there exist several other uncertainties, which 
may affect the extraction of the mixing matrix of scalar mesons. 

• Our results depend on the form factors calculated in the perturbative QCD approach, which is based on the fcy 
factorization. The application of PQCD approach to many channels involving the s-wave mesons is successful. 
Phenomenologically, the large annihilations can explain the correct branchin g ra tios and direct CP asymmetries 
of _B" 7r+7r~ and K~tt^ 31|, the polarization problem of i? ^ (pK* 32|, etc. Up to now, this approach 

is only at the leading order in l/ms expansion, without the next-to- leading order corrections. Since the scalar 
glueball is not very light, the power corrections proportional to msJniB may also be important. Moreover 



although this approach has been proved for the B ^ -k form factor |33j , there is no proof for the factorization 



of the i?-to-glueball form factor and this may introduce some potential uncertainties. 

• The hadronic inputs give another origin of the uncertainties. The decay constant of a scalar glueball may 
introduce sizable uncertainties (roughly 20%). Since there is not any theoretical study on the LCDA of scalar 
glueball, the higher but unknown Gegenbauer moments in the LCDA may also provide sizable contributions. 

• In the presence of mixing, the uncertainties from ordinary scalar mesons also affect the production rates of 
scalar glueballs sizably. For example in the two different scenarios for scalar mesons, the form factors for the 
quark states even have different signs, because of the different signs in the decay constants of ordinary scalar 
mesons 20|. Accordingly, the interference between quark content and giuon content differs in signs in the two 
scenarios. Moreover, the magnitudes of decay constants could also suffer from large ambiguities. We take the 
decay constant of /o(980) as an example. Based on the two-quark assumption for the /o(980) meson, the early 



study in the QCD sum rule [35| gives 

4(980) = (0.18±0.015)GeV, 



while the recent study |20[ predicts two sets of values 

//o(980) = (0.37±0.02)GeV[(0.46±0.025)GeV], 

at the scale p = IGeV [p = 2.1GeV). The decay constants for the isosinglet mesons near 1.5 GeV may also 
suffer from similar uncertainties and the interference between different components will be complicated. 

These quantities will inevitably affect the production rates of scalar mesons, and the extraction of the mixing ma- 
trix will require more precise, both experimental and theoretical, studies. For example, the experimental study on 
semileptonic B ao(1450) decays is helpful to constrain the contributions from the quark content. On the theoretical 
side uncertainties caused by the hadronic inputs could be reduced in a systematic and comprehensive study on the 
ordinary scalar meson and the glueball state. This is beyond the scope of this work and will be reported elsewhere. 
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D. Comparison with other works 

In Ref. the authors nwestigated the inclusive B XgG decay, where the large production rate of /o(1710) 
is viewed as the signal of a scalar glueball. The rate of B decays into a scalar glueball through gluonic penguin is 
expected to be sizable, but B — > XsG is an inclusive mode, which is more difficult than exclusive B meson decays 
experimentally. The background from the other mesons may also pollute their method. The authors in Ref. Q used 
symmetries of the penguin contributions and studied the production rates of a scalar glueball in charmless three-body 
B decays. Irrespective of the validity of neglecting contributions from tree operators (with large Wilson coefficients) 
and the uncertainties from the symmetry breaking effect, much more data on the B ^ S decays is required to make 
precise predictions. 

In Ref. 0, the author studied the exclusive channel B GK{K*) which could be experimentally detected. Since it 
is a nonlcptonic decay, the uncertainties would be larger than those in _B — > Glv and Bd,s Gl^l^ decays. B GK 
is a penguin-dominated process, which is purely induced by the loop effect in the standard model. The relevant 
Wilson coefficient is small and the result is sensitive to the next-to- leading order correction and/or the potential new 
physics effect. Hadronic uncertainties are also typically large in nonlcptonic decays. Contributions from the _B — > G 
transition form factor arc neglected in their analysis. The light-cone distribution amplitude is not used in their work, 
thus the internal structure of a scalar glueball could not completely reflected. It would be interesting to reanalyze 
this nonlcptonic decay channel by combining our prediction on the form factors and the other amplitudes studied in 
their work. 

Compared with the studies in the literature, we can see that albeit the production rates of a glueball state in 
semileptonic B decays arc not very large, the chosen channels B ^ Glv and Bg Gl^l~ are easily measured and 
rather clean. In inclusive decays B — > XsG decays, it is almost impossible to study the mixing between ordinary 
scalar mesons and a glueball state. 



E. Some discussions on nonleptonic B decays 

Semileptonic B decays are clean but in B — > f[)lv, the neutrino is identified as missing energy and the efficiency 
is limited; while the Bs — >■ fol^l~ has a small branching ratio. In these decays, the lepton pair does not carry any 
SU(3) flavor and the decay amplitudes receive less pollution from the strong interactions. The lepton pair can also 
be replaced by a charmonium state such as J/ip since J/?/' does not carry any light flavor either. B — >■ J/i/'/o decays 
may provide another ideal probe to detect the internal structure of the scalar mesons. In £? — > J/ipfo decay, the ss 
component will not contribute at the leading order in a^. For example, the B — >■ J/ipcj) decay has been set a very 



stringent upper limit [3J|: B{B — >• J/tpcj)) < 9.4 x 10"''. Thus B — > J/V'/o decay can ffiter out the glueball component 
and the fin component of a scalar meson. Meanwhile in Bg — > J /ipfQ decay, only the ss and the gluon component 
contribute. Moreover, the final state mesons in these channels are easy to reconstruct and these channels could have 
sizable branching fractions. If we use the factorization method, decay amplitudes are given as 



A{&^J/i^h) = -^VMi2e%,-PB)fj/^a2Ff^f"im%^). (35) 



Gt 

The Wilson coefficient 02 can be extracted from the B — > J/ipK decays [2f 

B{B° J/tpk^) = (8.71 ± 0.32) X lO""^. (36) 
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c 





(b) 



FIG. 3: Feynman diagrams of B — s> foD and Bs — > foD decays. These channels suffer large uncertainties. 



The branching ratios are roughly predicted as 

B{B" ^ J/^jfo{nn)) 



J (23^11) X 10-^ SI 
(1.5 ±0.5) X 10"^ 



(37) 
(38) 



where we have assumed the same dependence for all form factors and F^~^^ (0) = 0.3. The uncertainties are from 
the experimental data for B{B — > J/ipK) and the B ^ S form factors at the = point. For the Bg decays, the 
branching ratios are comparable with that of _B — > J/ipK: 



B{Bs -> J/4'Mss)) 
B{Bs J/i'G) 



r (6.5ttg) X 10-4 

\ {i.htll) X 10-3 

(2.4 ± 1.0) X 10-^ 



Scenario 1 
Scenario 2 



(39) 
(40) 



Such large branching fractions offer a great opportunity to probe structures of scalar mesons. With the available 
data in the future, the mixing problem between the scalar mesons will be solvable and the glueball component can be 
projected out in principle. 

If the power-suppressed annihilation diagrams are neglected, the charmful decays of B meson, B /qD, can also 
be used to constrain the mixing between scalar mesons. For instance in B^ — > Jq, the nn and gluon component 
contribute but the ss component does not, while in Bs — >■ Jq^ the nn component will not contribute, as shown in 
Fig. 131 Thus the mixing coefficients can also be determined if these two channels are experimentally measured. It is 
necessary to point out that this method may suffer from sizable uncertainties of annihilation diagrams [l8| . 



F. Glueball production in Be decays 

The ordinary light scalar meson is isospin singlet and/or flavor SU(3) singlet, while the glueball is flavor SU(6) 
singlet. Therefore it is difficult to distinguish them by the light u, d and s quark coupling. However, the light ordinary 
scalar meson has negligible cc component, while the glueball have the same coupling to cc as that to the uu, dd or ss. 
A clean way to identify a glueball is then through the cc coupling to the glueball. 

In B decays, the initial heavy meson contains a light quark, thus contributions of the gluon component always 
accompany with the quark content nn or ss. It is not easy to isolate the gluon content. The situation in the doubly- 
heavy B(. meson is different: it contains a heavy charm antiquark. The semileptonic Be — > foli^ decays would happen 
only through Fig. [IIa)(b) and (c) but not through Fig. [IJd) and (e). The observation of this decay channel in the 
experiments will surely establish the existence of a scalar glueball. Moreover the CKM matrix element in this channel 
is Vcb, thus the Be — >■ folv will have a sizable branching ratio. This channel will depend on the Be G transition 
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form factor which requires the less-constrained Be meson's hght-cone distribution amphtude. But even if we assume 
the form factor of i?c ^ G is smaher than the Be r]c form factor by one order, branching ratios of Be Glv{lv) 
are suppressed by two orders 

B{Be Glv) - 1% X 0.01 = 10-^ (41) 

where the branching ratio of Be — > rjclv has been taken as 1%. Although more quantitative results need the precise 
Be ^ G form factors, which requires the Be wave function with large uncertainties, the order of magnitude of branching 
ratio is convincing. This branching ratio is large enough for the experiments. One only needs to reconstruct the /g 
scalar meson in the final state and also the Be meson mass in the intermediate state, so that to make sure that 
the scalar meson is produced from two gluons. That experiment is achievable even if the /g meson is not a pure 
glueball, but at least has a large portion of it. The disadvantage here is the missing neutrino in detector is hard to be 
reconstructed in hadronic machines, like LHCb. A future Z factory is an ideal place for this channel to be measured 

Be — > foT^~ is another potential mode to figure out the gluon content. But in this mode, the fin component also 
contributes through the annihilation diagrams. The b and c quark annihilates and the d and u quark are created. The 
CKM matrix element Veb and the Wilson coefficient ai are the same with the emission diagram for the i?c-to-glueball 
transition. The offshellnes of the two internal particles in annihilation diagrams are of the order . The electrowcak 
vertex is the V ~ A type and the decay amplitude is proportional to the light quark mass. Thus the decay amplitudes 
via annihilation diagram for the hn component are expected to be suppressed. As a result, the Be — > /o^r^ also filters 
out the gluon component of the scalar meson as an approximation. 

V. CONCLUSIONS 

In this work, we have investigated the transition form factors of B meson decays into a scalar glueball state in 
the PQCD approach. Compared with the form factors for the quark content, we find that form factors for the gluon 
content have the same power counting in the expansion of Aqcd/'tis but they are numerically a little smaller. The 
pure glueball states can be detected in semileptonic B Glv and Bg Gl^l^ in the future, since they have a 
sizeable branching ratio. 

If a scalar meson is a mixture of a glueball and an ordinary meson, we investigate the possibility to extract the mixing 
mechanism from semileptonic B decays, such as S — >■ f[)lT> and Bg — )- fol^l~ decays. The nonleptonic B — > J/ipfo 
and Bg — > J/^pfo decays are also analyzed. To be specific, we discussed the production rates of scalar mesons under 
two different mixing mechanisms, and we find that the differences in B and Bs decays are helpful to distinguish the 
two mixing mechanisms. To avoid the interference between the quark and the gluon component, we find that the 
Be ~> fali^ and Be — > /ot" will project out the gluon component of a scalar meson cleanly. Our results can be 
generalized to the other glueballs. 
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Appendix A: hard kernels 



Hard scales and functions in the B ^ G form factors are defined by: 



ta 

tb 

tc 

Ea{ta)ha 
Ec{tc)hc 

where Zi is defined as 



max [y/pxiTUB, 1/61, 1/&2] , 

max [^/pxirriB, y/px^ruB, I/&1] , 

max [\/p?7iB, 1/61, 1/62] , 

+e{h2 - bi)Io{y^px^mBbi)Ko{y^px^mBb2)] Ko{y/px^mBb2), 
as{tb)c-x:p[-SB{tb) - Scitb)] I dz 

Jo ^V'^l 

{Ea{ta)ha)\ 



Ki{^Z^bi) 



Zi = px2m%z + p.T2.Tim|(l - z) 



The Sudakov factors are given by: 

SB{t) 

Scit) 



rriB , 
s I xi—^,bi 



dp 

l/bi M 



V2 



So 1 X2 —,b2 ) + SG (1 - X2) —,b2 



V2 



where the exphcit form for the function sg{Q, b) is 

rQ dp 



SGiQ,b) 



l/b M 



hi 



A{a{p)) 



A = — Ca, 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 

(A6) 
(A7) 

(A8) 

(A9) 
(AlO) 

(All) 



with C'a = 3. Since there is no endpoint singularity in these diagrams, the contribution from the endpoint region 
will not be large and we will neglect the jet function arising from the threshold resummation. We adopt the one-loop 
expression of the running coupling constant a^, when evaluating the above Sudakov factors 
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